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ABSTRACT 



In this Letter, we propose that a microphysical process takes a vital role 
in the shocked region in which the prompt emission of GRBs is emitted. The 
turbulent energy is included in the internal energy transferred by the kinetic 
energy of the shock. It dissipates through stochastic acceleration for the electrons 
to supply the early X-ray emission in the phase of shallow decay. We put the 
constraints on the time evolution of microphysical parameters. The early X-ray 
fluxes can be obtained by this scenario and these results are consistent with the 
Swift observation. 

Subject headings: gamma rays: bursts — radiation mechanisms: general — X- 
rays: general — turbulence — stochastic acceleration 



1. Introduction 



In the Swift era (IGehrels et al.ll2004j ). the X-ray telescope (XRT) observation (IBurrows et al. 



20051 ) provides the complete light curves of gamma-ray bursts (GRBs) in the 0.2-10 keV 
band. One of the most interesting discoveries is the so-called shallow decay segment: 
the flux plateau of F cx t~°[ 5 within 10 3 — 10 4 second after the trigger (ICampana et al. 



2005 



De Pasauale et a 



Liang. Zhang fc Zhang 



2006). Recent statistic analyses (INousek et al.ll2006l ; lO'Brien 



2006 



200 



?]) have revealed that the phase of shallow decay in the X-ray af- 



terglow might be a common feature of the long GRBs. The shallow decay in the early X-ray 
light curve is still a mystery, although theoretical explanations have been put forward from 
several aspects (see Zhang 2007 for a compr ehensive review). Most of the models are: hydro - 
dynamics of the sho ck by energy injection (e.g.. Granot &: Kumarll2006l ; IZhang et al.ll2006l ) 



geometry of the jet Je.g.. Eich 



parameters ( Fan fc Piranl 120061 ; iGranot. Konigl fc Piran 



er h Granot 



2006; 



Toma et al 



2006 



2006h~ varying microphysical 



Panaitescu 2006: Ioka et al. 
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2006 1) , late promp t emission (jGhisellini et al.l 120071 ) or up-scattered forward-shock emission 
f|Panaitescull2007h . 



From the point of microphysics in the hydrodynamic ev olution, the nature of co upling 
between the electrons, protons and magnetic field is complex (j Chiang fc Dermerll 19991 ). Usu- 
ally t he simple way is to assume the equipartition between electrons, protons and magnetic 
field ( iPanaitescu fc Meszazrosll2006l ). However, in this Letter, we propose that the kinetic 
energy of the relativistic shocks in the plasma has been converted to the internal energy as 
three parts: (1) eb, which is the energy of magnetic field; (2) e e and e p , it means that the 
electrons and protons/positrons are accelerated by the shocks, normally, this is the process 
of first-order Fermi acceleration; (3) e t , which presents the turbulent energy, and this energy 
would sustain a relatively long time (see Section 2). The last part has not been taken into 
account by the former research. In our novel scenario, the relativistic electrons accelerated 
by the first-order Fermi acceleration emit the gamma-ray by synchrotron radiation, after the 
decrease of the prompt emission tail which is shown as the deep decay in the early X-ray 
light curve, indicating that the internal shocks are abated, the follow-up turbulence and its 
effects might be dominated. The turbulence could transfer its energy to the electrons via 
second-order Fermi acceleration, which is also called as the stochastic acceleration. There- 
fore, in this turbulent region, due to the resonant interaction between electrons and plasma 
waves, the electrons buried in the magnetic field are re-accelerated by the stochastic accelera- 
tion. The emission of these electrons dominates the shallow decay phase, until the turbulent 
energy dissipates and the external shock sweeps the surround medium thus the deep decay 
appears again. 

In Section 2, we review the Fokker-Planck equation and list the coefficients associated 
with the turbulent term. In Section 3, the turbulent parameter e t , as same as e e and eb, 
is introduced. Due to the turbulence, these microphysical parameters evolved with time are 
constrained by the process of stochastic acceleration. Finally, we select these relations to 
reproduce the feature of shallow decay. The discussions are given in Section 4. 



2. Stochastic Acceleration 



The stochastic acceleration was suggested as an non-neglected mecha nism to produce 
the high energy particles in GRBs (IWaxmanl Il995l ; iDermer fc Humill200ll ). The numerical 
simulation has confirmed that the stochastic acceleration in the relativistic shocks plays a n 
important role on the particle energy distribution and evolution (IVirtanen &: Vainioll2005l ). 



In general, the charge particles are expected to be accelerated through resonant interac- 
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tions with the magnetized plasma waves. For this stochastic acceleration, particle distribu- 
tion function /( 7 , t) = dN(j, t)/dj satisfies the Fokker-Planck equation, which is the kinetic 
equation of single electron in the energy space: 



0/(7,*) 
dt 



^[U( 7 )/( 7) t)]_ A[A( 7 )+J5 L /( 7 ,t)] 



T{l) 



+ Q(i,t) 



(i) 



where A (7) is the acceleration rate, .0(7) the diffusion rate, El energy loss rate, T(j) the 
mean escape time of electrons and Q(7, t) the source term. The turb ulent spectrum is usually 
empl oyed by the form of Kolmogorov or Kraichnan as W(k) oc k~ q (je.g.. Zhou fc Matthaeus 
1990h . 



Park and Peterosianl (119951 ) already explored the various analytic solutions and they 



illustrated that the coefficients of the equation are related to the microphysics of turbu- 
lent plasma. In that paper, they offered one special example, the so-called hard-sphere 
approximation, w here the index of turb u lent s pectrum q = 2. This detailed work has been 
complemented by iBecker. Le fc Dermerl (120061 ). in which the time-dependent Green's func- 
tion was reexamined and the case of q < 2 was considered. However, the energy loss, such 
as synchrotron or inverse Compton, was not included in the solution. 



Wang fc Mad (120011) calcu l ated t he equation which is applied in the spectral variability 
of blazars. Recently, iLiu et al.l (120061 ) studied the stochastic acceleration in Sagittarius A* 
adopting q = 2. Under the physical cond it ions of Sagittarius A*, the stead y state of f(j) 
was derived (ILiu. Melia fc Petrosianl 120061 ). iManolakou. Horns fc Kirk! (120071 ) have obtained 
the time dependent solution of /(7) taking into account the cooling via synchrotron and 
inverse Compton radiation. In this Letter, the calculation is simply on the coefficients and 
we avoid to give any details about fipf). 

In order to explain the shallow decay phase of GRB X-ray afterglow, assuming (i) the 
scatter path is much less then the length of the turbulent region, (ii) there is no continual 
injection and the accelerated electrons have the synchrotron cooling thus the steady state 
exists, we have the simplified Fokker-Planck equation of df /dt = and all the coefficients 



are time-indepen dent . Here, we repeat the formulae produced by lDermer. Miller fc Lil (119961 ) 
and IWaneJ J2002h : 

^-(Ub + U^P) 2 , (2) 



E, 



Mi) 



3m e c 
irq-lU t 



2 q U 



fickirLky-y- 1 , 



B 



~ 1 Ut -p 2 g ck(r L ky- 2 p^-\ 



2q(q + 2)U B ' 



(3) 
(4) 
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T(j) = *- q —^ck{r L kY-^-\Mfp-\ (5) 
z q Ub 

where (3 = v/c, 7 is the Lorentz factor of electron, k = (ct) _1 the wavenumber, tl = 
m e c 2 /eB the nonrelativistic Larmor radius of the electron, ot the Thomson cross section, 
p dimensionless momentum. Ub, U t and U p h are the magnetic energy density, turbulent 
energy density, and photon field energy density respectively. Suppose the the particles are 
accelerated by Alfven turbulence, the group speed (3 g is equal to Alfven speed. 

In the traditional models, two microphysical parameters, the ratio of magnetic field 
energy to the total internal energy Eb = Ub/U and the fraction of the total internal energy 
that goes into the random motions of electrons e e = U e /U, are drawn into the study of GRBs 



( le.g.. Sari. Narayan fc Piranlll996l ). In our scenario, the turbulent energy is included in the 
total internal energy. Similar to the e e and Eb, we define the dimensionless parameter e t to 
present the turbulent property as: 

et = jj. (6) 

Since the dominant emission is synchrotron radiation during the shallow decay phase, 
the term of photon field U p h which is associated with the inverse Compton process can be 
ignored. Of the single electron, the stochastic acceleration rate is balanced to the energy 
loss rate through synchrotron radiation: 

A{l) = El- (7) 

The energy ratio of accelerated electrons in the turbulent region is equal to the total 
stochastic acceleration rate: 

dee = [A^fj^dq 

dt j7( 7 )d 7 • 1 > 

The dissipation timescale of turbulent energy can be estimated by: 

m EtU alls , s 

< E L > dE e /dt 

where a = Et/Es- The relativistic electrons re-accelerated by stochastic acceleration are 
assumed to have energy distribution f(j) oc 7~ p with the limits of [jmin, Jmax] and produce 
X-ray emission. For p = 2.2, the timescale Td is given by: 

T d = 8.6 x 10 3 x a x ( -L)-^(^_)°-2( " )0.2 ( 7mm ) - 1 . 2 (1Q) 

v ioo ; v o.oor v io 2; v io.(r ' v ; 

where T is the bulk lorentz factor, 7 m j„ is the minimum lorentz factor of the accelerated 
electron by turbulence and n is the number density of shocked medium. The dissipation 
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timescale of the turbulent energy is coincident with the observational timescale of shallow 
decay phase. 

Therefore, in the shocked region with turbulence, the evolution of microphysical param- 
eters e e , Eb and e t can be derived from the coefficients of Fokker-Planck equation. From 
equation (7) and (8), we obtain: 

^oc^ 1 ^ (11) 

and 

1 oc ^g-^V"*. (12) 



dt 

Furthermore, in this turbulent region, we assume that these microphysical parameters 
have the forms of e e oc t a , eb oc t b and e t oc t d . We insert them into equation (11) and (12) 
and put the constraints on a, b and d by two algebra equations: 

b = 2d/q + 2(1 -q)/q (13) 

and 

a-l = 6(2-g)/2 + (l-g) + d. (14) 

We choose the turbulent spectrum of Kraichnan with the index q = 4/3. All of the possible 
values of a, b and d are shown in Figure 1, 2 and 3. 

From the standard afterg low model, the fluxes in the early X-ray band were written by 



Sari. Piran fc Narayanl (119981 ) . There are two different limits: adiabatic and radiative case. 



For explanation of the emission in shallow decay phase, we obtain the early X-ray light curve 
under the adiabatic case as: 

F OC e P- 1 £: (P- 2 )/ 4 £( 2 - 3 P)/ 4 oc t (b P -2b-3p+2)/4+a(p-l)^ ^ 

As an example, we adopt the index q = 4/3 and p = 2.2. Four evolutionary fluxes are 
achieved in the table [TJ they are corresponding to the minimum and maximum values of a 
and b. While for the radiative case, the X-ray flux is: 

F OC eP-l e CP- 2 )/ 4 ^(2-6p)A oc £&(p-2)/4+(2-6p)/7+o(p-l) _ 

Table [2] lists the possible fluxes evolved with the time. 



3. Discussions 



Due to the turbulent energy and its dissipation in the shocked region, after the gamma- 
ray radiation, the electrons are re-accelerated by the stochastic acceleration and emit the 
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early X-ray fluxes. Therefore, the energy injection in the central engine is not required to 
reproduce the emission in the shallow decay phase. Our results manifest that the temporal 
feature in the shallow decay phase represents the process of turbulent dissipation. The 
microphysical parameters, e e , eb and e t , are varied with the time in the shocked region. 
Given the special values of p = 2.0 and a = 0.5 in the adiabatic case, our calculation can 
represent the typical observational flux as F oc t -0 ' 5 . Moreover, it is noted that the same 



shallow decay phase is also detected in the optical band as well (je.g.. Mason et al.l 120061 ) . 
our simple interpretation is that the synchrotron emission of X-ray band and optical band 
may be original from the same shocked region wit h turbulent energy. Another advantage of 



our model is that the crisis of radiative effici ency (ILloyd-Ronning &: Zhang! 12004 ; lloka et al. 



20061 ; iFan fc Piranl 120061 ; IZhang et al.l 120071 ) is therefore dispelled without any additional 
assumption of ejection/ejecta from the central engine. From this point of view, we support 
the internal shock pattern of standard fireball model. 

In this work, we assume /(7) oc 7~ p where p = 2.2 as the typical value. The detailed 
calculation of the energy distribution f(j) is not needed, because in our scenario the timescale 
of acceleration is much smaller than that of turbulent energy dissipation. Since the timescale 
of acceleration is also smaller than that of shock hydrody namics, therefore, during the p hase 



of shallow decay, the index of spectrum does not change (ILiang. Zhang fc Zhang 
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Generally, the varied values of e e a nd Eb lead to the diff erent hydrodynamic evolution 
and emissions in the entire afterglow ( Mao fc Wang 2001al lbh . In this Letter, the early 
X-ray fluxes in the shallow decay phase are estimated either in the adiabatic or radiative 
case through the turbulent process. The average slope of radiative case is deeper than that 
of adiabatic case. Compared with the XRT observation, the hydrodynamic evolution of 
adiabatic case could be the realized regime during the shallow decay phase. Although it is 
hard to obtain the values of e e and Eb in our scenario, we put the constraints of e e <C 1 and 
£b<1. More observational samples are particularly requested for the further investigations. 
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Fig. 1. — The a-b diagram. The related values of a and b are shown in the region noted by 
the capital P. 
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Fig. 3. — The b-d diagram. The related values of b and d are shown in the region noted by 
the capital P. 
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Table 1: The limit fluxes of shallow decay for the adiabatic case 



adiabatic 


case 




Q>max 0.5 




-2 


F oc t" 1 - 25 


F cx t" - 65 


b-max — 


-0.5 


F oc r 1 - 175 


f oc r - 575 



Table 2: The limit fluxes of shallow decay for the radiative case 



radiative case 


O'min 


Q>max 0.5 


b ■ — —2 

u min ^ 


f oc r 1 - 5 


F oc r 09 


h - — n ^ 

u max u ■ u 


f oc r 1 - 575 


F oc r - 975 



